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Required Background 



Required Background 

1. Thermodynamics (PHYS253) 

(Covers most of the basic ideas. It would help to look 
through the notes again.) 

2. Quantum and Atomic Physics (PHYS255) 

(You need: Schrodinger's equation, particle in a box, simple 
harmonic oscillator, Zeeman effect, fermions and bosons.) 



Basic Statistical Mechanics 



2 



Recommended reading 



1. Statistical Mechanics - A Survival Guide, 
A. M. Glazer and J. S. Wark 

Oxford University Press, 2001 

2. Basic Superfluids 
Tony Guenault 

Taylor &l Francis Inc., 2003 

Both available as ebook in the Liverpool University library. 
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Thermodynamics, briefly. 
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A brief review 



1. Thermodynamics 

2. Quantum mechanics 
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4. Paramagnetic salts 

5. Ideal gas 
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Heat capacities 

This graph shows the heat capacity of copper. At high 
temperatures (T), it tends to a constant. At low T, it falls to 
zero. 




The main theme in Part 1 is to learn how to explain such 
behaviours. 
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Entropy 



The study of heat engines in the 19th century led Rudolf 
Clausius to propose in 1982 the second law of thermodynamics 
in this form: "In a closed cycle, the sum of 5Q/T is zero." 

http : //en . wikipedia . org/wiki/Entropy 

5Q is an infinetesimal heat change. Later, 5Q/T came to be 
written as 5S, or 

dQ = TdS, 

and S became known as the entropy. 

Clausius statement implies that this abstract quantity S is a 
state function. E.g. if you compress and expand a gas and go 
back to the initial p and V at the end, then S also returns to 
the same value. (Note that Q itself is not a state function.) 
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Quantum mechanics, briefly. 
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Chemical potential. 



The chemical potential fi is equal to the Gibbs free energy at 
equilibrium. The Gibbs free energy is 

G = U + P V - TS. 

One simple use of n is to find the energy change between two 
phases at constant p and T: 

Afj, = /\U + pAV - TAS. 

At phase transition between liquid and gas, this would be the 
energy difference between 1 mole of liquid changing to gas, and 
1 mole of gas changing to liquid at the same T and p. So 

A^i = 0 

gives the condition for phase equilibrium, because then there is 
no energy advantage in more liquid changing to gas, or the 
other way round. 

http : //en . wikipedia . org/wiki/Gibbs_f ree_energy 
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Schrodinger's equation. 



This is the ID Schrodinger's equation. 

^ + V( X )i, = ETp, 
2m ax 

where m is the mass of a particle, ip is the wavefauction, V(x) 
the potential it experiences, and E its energy. 

In free space, V(x) would be zero. Solving the equation then 
gives 

^ _ e ikx Qr e -tkx^ 

where k is related to E given by 

2m 
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Momentum. 



The Scrodinger's equation is based on the postulate for the 
momentum p x of a particle: 

-in— = p x - 
ax 

We use this on a wavefunction ifi like this: 

-in— = p x xp. 
ax 

Applying this to the free space wavefunction e tkx , we find 

p x = hk. 
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Statistical mechanics 
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Particle in a box 



In a ID box of length L, the potential V(x) could be zero for x 
between 0 and L, and infinity otherwise. 

This means that the wavefunction must be zero outside the 
box, and also zero at the "walls" at x equals 0 and L. 

Inside the box, it is like free space, so t/> is e lkx or e~ %kx . In order 
for ip to be zero at the walls, these can be combined to give 

ip = sin kx = . 

2i 

Applying the condition that this is zero at the walls, we find 

nir 

k = T> 

where n is an integer. This shows that the energy is quantised: 

_ h 2 k 2 _ n 2 Tv 2 h 2 
2m 2mL 2 

Basic Statistical Mechanics 13 



Microstates 



In quantum mechanics, energies are often quantised. We have 
seen that this happens for a particle in a box. 

Consider a system with N particles. This could be atoms in a 
solid. For now, think of each one like a particle in a box. 

Suppose every particle has energy levels ex, £2. ■■■ 

A microstate refers to a specific arrangement of particles 
among the energy levels, e.g. particle 1 in e 3 , particle 2 in e±, ... 

A macrostate refers to a distribution of particles among the 
levels, e.g. 6 particles are in e\, 4 particles are in e 2 i •■■ 



Basic Statistical Mechanics 



15 



Macrostate. 



Suppose there are n\ particles in s\, n 2 particles in e 2 , ... 

The macrostate can be written as (n 1 ,n 2 ,...), or simply 
where i = 1, 2, ... Then 

AT = ni + ri2 + ••• 

Suppose the system is thermally isolated from the surrounding, 
so that the total energy E is constant. E is given by 

E = n\e\ + n.2£2 + ••• 

These two equations impose constraints on the number of 
particles possible on each energy level. 
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Constraints 

Assuming the constraints 

N = n\ + ri2 + ••• 

and 

U = n\e\ + n2£2 + ••■ 

it is possible to solve for the most probably macrostate using 
the method of Lagrange multipliers. 

Before that, we need to write down the number of possible 
microstates for a the macrostate. This is just the number of 
ways to arrange N different objects in some boxes, with n\ 
objects in box (level) 1, n 2 objects in box 2, and so on. 

Using the method of combinations, the number of 
arrangements is 

n = ^_ 

n\\n2\... 

The problem is to find the values of n, that gives the largest Q. 
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Postulate of equal a priori probabilities 



More than one microstate can have the same macrostate. E.g. 
if the macrostate has 4 particles are in e 2 , then these 4 
particles can be any of the N particles. 

Basic postulate of statistical mechanics: 

"A system is equally likely to be in any microstate that satisfies 
the constraints." 

This means that a system is most likely be be in the 
macrostate with the largest number of microstates. We may 
call this the most probably macrostate. 

Using the method of statistics, it can be shown that if N is 
very large, then the probability of other macrostates is very 
close to zero. (E.g. when you take the average of more 
readings in an experiement, the error gets smaller.) 
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Method of Lagrange multipliers. 
The steps are as follows. 

1. Write down the Lagrange function In Q + X^N + X 2 U. 
Two new variables, X 1 and A 2 , have to be introduced. 

2. Differentiate with respect to each of and A 1( A 2 . 

3. Set all derivatives to zero. Solve the equations for n^. 
The answer is: 

ni = Aexp(X 2 Si) 

where A is exp^). 

There are two unknowns A and A 2 . In principle, we can 
substitute into the above constraints and solve for these. 
Unfortunately, the equations are too complicated. 
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Boltzmann's postulate 



Around 1872, Ludwig Boltzmann postulated the connection 
between entropy S and the number of macrostates Q: 

S = k B \nn. 

This formula allows us to use thermodynamics, e.g. using 

dQ = TdS, 

to find the answer to the most probable macrostate. Assuming 
Boltzmann's postulate, it can be shown that 

A 2 = -l/k B T, 
so that the distribution becomes 

rii = Aexp(-ei/k B T). 
This is now called the Boltzmann distribution. 
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Paramagnetic salts. 

Heat capacities tend to increase with temperature, as particles 
can reach higher energy levels. 

At very low temperatures, paramagnetic magnetic shows a 
peak in its heat capacity. 




Tamperature 77 K 

J.W. Stout, W.B. Hadley, J. Chem. Phys. 40, 55 (1964) 

This unusual behaviour is called the Schottky anomaly. 
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Paramagnetic salts 
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Schottky anomaly 



The Schottky anomaly happens for paramagnetic materials in a 
magnetic field. At very low temperatures, there is little 
contribution from electrons and phonons (lattice vibration). 

Consider a material with ions that have spin 1/2. This is called 
a spin 1/2 paragmagnet. When a magnetic field B is applied, 
the energy level of each ion split into two. 



0 

s, = -fS 



e 2 =pB 



it: 



Ae = 2fiB 



The energy levels are 

e = -fi.B = -/j-bB or + MB-B 
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Lets see if we can "deduce" the Schottky anomally. At 0 K, all 
electrons are at the ground state, so the total energy U is 
-Nfi B B, where N is the number of the ions. 



At high temperatures, the thermal energy the electrons is much 
larger than the difference between the two magnetic energy 
levels. Then the electrons are equally likely to be in -fi B B or 
+HbB, so U is zero 

0 T 
U 

-Nu B B 

A sketch of U against T would look like this. 
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Now that we know the gradient is zero at both 0 K and high T 

o L_ 

u 

-Nn B B 

we can sketch the heat capacity, which is the gradient of U. 




T 

We have just deduced the Schottky anomaly. 
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Heat capacity 



In principle, we can now find the heat capacity C using 



However, the gradient of U neat 0 K is not known. We can find 
this using the Boltzmann distribution. 

Near 0 K, nearly all electrons at at the lower state. The very 
small fraction in the higher state is approximately given by the 
Boltzmann factor exp(-e/fc B T) (why?), where 

e = 2 Mj5 B. 

The total energy is then 

U = -N^ B + Nfj. B exp(-e/fc B T). 

Differentiating this shows that the gradient goes to zero as T 
goes to 0 K. 
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We shall now use the Boltzmann distribution to derive a 
formula for U and C. There are two levels: 

e\ = ~hbB and £2 = +/ i _B^3- 
The total energy is 

U = J2 n i £ i 

The number of electrons in level i is given by 

m = Aexp(-£i/k B T). 
We must solve for A first. The total number is 

N = A£exp(- £l /fc B T). 
The sum is called the single particle partition function: 

Z SP = J2 ex P(- £ i/ k BT) 

Solving for A gives 
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So 

N 

rii = —— exp(~ei/k B T). 
Z SP 

The total energy is 

U = Yl n i £ i 

Substituting nf 

N 
Z SP 

Substituting for ef 

N 

U = —— [-n B Bexp(n B B/k B T) + n B Bexp(-n B B/k B T)]. 

Z SP 

and in Z$p as well: 

_ n -»bB exp{^i B B/k B T) + [i B B exp(-fj, B B/k B T) 
exp(fi B B/k B T) + exp(-v B B/k B T) 
Simplifying, we find U: 



U = -N^ B Btanh(fx B B/k B T). 
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Differentiating U gives the heat capacity: 

C = Nk B f^Vsech 2 (^ B ). 
which also gives the Schottky anomaly that we deduced: 
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The formula 

U = ~Nfj, B Btanh(fi B B/k B T). 
does indeed give the graph that we deduced earlier on: 



-Nu B B 




We note here that a useful formula from the above formula for 
rij, is the probability that an electron is at level ef 

= iH = exp(sj/k B T) 
Pl N Z 
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The entropy. 



We can find the entropy using 



dS = 



dQ dU 



T T 

In our case, Q is just U since no work is done. So 

• dU r CdT 

This can be done analytically using the earlier formula for C. 
The result is 

f li B B^ 



= /f = / : 



S = Nk B \n 



2 cosh 



k B T 



T \k B T ) 
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The entropy. 



The graph for the entropy looks like this: 

S 



Nk B ln 2 




Again, we would be interested in the high and low temperature 
behaviours. 

Recall the formula for entropy: S = k B In Q. 
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Ideal Gas 



Basic Statistical Mechanics 



34 



Behaviour of the entropy 



At low temperature, most atoms would be in the ground state. 

+H B B 

-H B B — tf— 

There is only 1 way to arrange this. So the entropy would tend 
to S = k B \n 1 = 0. 

At high temperature, the difference between the energy levels 
would be small compared with the k B T in exp(— e/k B T). 

+n B B — ( — 

-H B B —J 

The atom is equally likely to be in either level. There are 2 
possible arrangements for each of the N atoms - 2 N in total. 

So the entropy would tend to S = k B In 2 N = Nk B In 2. 
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A particle in a 3-D box 

In an ideal gas, we assume that there is little interaction 
between atoms. 

From the kinetic theory of gases, we know that the average 
energy of each atom in an ideal gas is 

1 2 3, m 
—mv = —kaT. 

2 2 " 

This means that the total energy is 

U = -Nk R T 
2 u 

and the heat capacity is 

c = *L = 3 Nk B . 

dT 2 a 

http : //en. wikipedia. org/wiki/Kinetic_theory 

We shall rederive this here using statistical mechanics. The 
ideas will also be used for electrons and phonons. 
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Consider a particle in a 3-D box. 





Solving the Schrodinger's equation, we get the same 
quantisation condition as in the ID box. For the x direction, 

n x iv 



Likewise for the y and the z. 
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Density of states. 



This suggests that we may approximate the graph to a curve, 
and the sum to an integral. 




Within a small interval de, rij is nearly constant. Then the sum 
of the energies of particles is approximately n,ej times the 
number of states in this interval. 

If we divide this number of states by de, we can think of the 
the answer as a kind of density. This density of states (DOS) is 
often denoted by the function gfe). 



The number states in de is then g(e)de, and the energy in the 
interval given by nj£jg(£j)cte. 
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Density of states. 



To find the heat capacity of N particles, we need the energy 

U = J2 n i £ i- 

For the spin 1/2 salt, there are only 2 levels to sum. For the 
particle in the box, there are many levels. 

At room temperature, many of these levels could be occupied. 
This is because for an ideal gas, spacing between energy levels 
is much smaller than the average energy of a particle. (Prove 
it.) 

As a result, a plot of n, against e i may look like this. 
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Density of states. 



The total energy 

U = Y. n i £ i- 

can be written as 

Approximating this as an integral gives 

U = J n(e)eg(e)de. 
where we have replaced the subscript with the argument e. 

This creates a new problem. We must now find the DOS. 
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Counting states 



"he wavevectors 



kx 



a 



are quantised at uniform intervals in all three directions. 

we can imagine a fc space in which the the x coordinate is k x , 
and so on. If we use a point to represent each state, we would 
get a lattice like this. 
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DOS 



We find that the total number of states below fc is 

Vk 3 



G(fc) 



6tT 2 



where V is the volume a 3 of the box. 

Suppose k is increased by dk. Then the number of states 
increases by dG. So the density of states is 

dG(fc) Vk 2 

9 ( k > = av. = o~~2 ' 
dk 2tt z 

We have earlier defined the DOS in terms of energy, g(e). 



Important: To convert the variable to energy, it is not a simple 
matter of substituting the formula 

h 2 k 2 

e = . 

2m 
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Counting states 




First, we find the number of all states below energy e. This 
means the number of points within 1/8 (why?) of the sphere 
of radius fc, where fc is the wavevector for e. 

Since the spacing between points is n/a, the volume associated 
with each point is (n/a) 3 . Therefore, the total number of 
states is: 
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Changing variables 



The reason is that for the same interval e, the size of the 
interval dk can be different at different energy because the 
relation 

h 2 k 2 

e = 

2m 

is not linear. (Prove it.) 

The correct way is to use the method for probability density 
function: 

ge{e)de = g k (k)dk. 
http : //en . wikipedia. org/wiki/Probabi lit y_density_f unction 

The subscripts are added here to emphasise that g(e) and g(fc) 
are different functions. 
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Changing variables 



Rearranging, 



as 



We can now substitute the relation 

h 2 k 2 



2m 

to find ffe(e). The answer is 

9(e) = — ^( 2 ™ £ ) 1 

where we have dropped the subscript again, as is the normal 
practice. (Misleading, but you can tell from the argument.) 



For comparison: 



2n 2 ' 
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Ideal gas energy states 



In an ideal gas at room temperature, it is very unlikely for two 
or more atoms to occupy the same energy state. So we can 
assume that each state is occupied by only one atom. This 
simplifies the calculation of the macrostate. 

First, we should justify the assumption. 

At a temperature T, we know from simple kinetic theory of 
ideal gas that the average energy of a gas atom is about 
3k B T/2. 

We need to know how many energy states there are below this 
energy at room temperature. 
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Changing variables 



It is also useful to have the formula for total number states 
in terms of energy. 

Since no e or k intervals are involved, we can simply substitute 

h 2 k 2 



2m 

The answer is 
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Ideal gas energy states 



We can use the formula we have derived: 

GW = ^(2-) 3/2 
by setting e = 3k B T/2. 

Putting in the numbers for one mole of ideal gas at room 
temperature, we find that the number of states below 3k B T/2 
is about 10 30 . 

1 mole of gas contains about 10 24 atoms. 

This means we have about 10 30 4- 10 24 = 10 6 states per atom. 
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Ideal gas energy states 



This means we have about 1 million energy states for every 
atom. 



So it is extremely unlikely that two atoms would ever occupy 
the same energy state. 

We can use this result to derive the energy distribution of these 
atoms, with the help of the Lagrange multipliers. 
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Macrostate. 
Use the method of combinations, we find 

possible arrangements. 

To find the total number of arrangements for all intervals, we 
must multiply together the answer for every interval: 



Since the number of energy states g^ is much larger than the 
number of atoms n it we can simplify this to 
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Indistinguishable particles. 



Consider a interval Si around the energy e^. Suppose there are 
are atoms and g^ states in this interval. 

We need to find the number of microstates. We cannot do it 
the same way as in the paramagnetic salt example. This is 
because the particles in the ideal gas are indistinguishable. 

(Although the ions in the salt are identical, they are 
distingushable - because of their fixed positions in the solid.) 

The problem then is the same as that of arranging rn identical 
balls and g^ identical partitions in a row. 
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Most probable macrostate 
We can now apply the Lagrange multiplier method. 

The answer is 

rii = A gi exp(-£i/k B T), 
which is the Boltzmann distribution again. 

We can find A using 

To approximate to an integral as before, recall that the number 
of states gi in the energy interval is given by g(e)d(e). 

Then m above becomes 

n(e)de = Ag(e) exp(-e/k B T)de, 

and 

N = a[ g(e)exp(-ei/k B T)de. 
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Energy distribution. 



The integral 



Z S P = J 9(e) exp(-e/k B T)de. 



is just continuous energy version of the single particle partition 
function that we have used earlier (for the spin 1/2 salts): 

%p = E ex p(- £ */ fc s T )- 

So we get the same relation 

N 

N = AZ SP , or A = 



Z SP 

Substituting the formula for g(e), we find 

'27rmfc B T\ 3 / 2 



Z SP = V 



V \2nmk B Tj 



3/2 
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At the Columbia University in New York, 



Miller and Kusch measured the speeds of the atoms in 
potassium vapour. They agree very well with the formula. 




Fig, 2. Schematic diagram of the apparatus designed 
to measure velocity distributions. 




REDUCED VELOCITY 



Fig. 4. Typical potassium velocity distributions. The vapor 
pressures in tiie n vena are given in Table IT. Run 31 was made 
with thick oven slits and runs 57 and &0 with thin slits. 



This was published in Physical Review in 1955. 
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Maxwell-Boltzmann distribution 



and 



The energy is just the kinetic energy, e = rnv 2 /2. We can 
change to variable v using 

n(v)dv = n(e)de. 

This gives the velocity distribution 

/2\ 1 /2 / m \ 3 / 2 
n(v) = N \-J — exp(-mv 2 /2k B T)v 2 

This is the Maxwell-Boltzmann distribution. 
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When would it break down? 

We have assumed that: "It is extremely unlikely that two 
atoms would ever occupy the same energy state." 



Recall that we have arrived at this by looking at the number of 
states below the mean energy 3k B T/2, when T is at room 
temperature. 

If T is very small, there may be far fewer states below the mean 
energy. Then it would be very likely for two atoms to occupy 
the same state, and the Maxwell-Boltzmann distribution would 
not be valid. 

Instead, quantum statistics - such as Fermi-Dirac or 
Bose-Einstein - have to be used. We shall learn about this. 
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Partition function 



That there is a set of formulae that we can use to find the 
energy from the partition function. 

The energy is given by: 

U = k^-^ 

where the many particle partition function Z N for for 
indistinguishable particles is 

Z N 

' J N 



7 _ ^SP 



AT! 

(Prove it.) 

Substituting this, we find 



TT 3Nk B T . „ dU 3^, 

U = — and C = — = -Nka, 

2 dT 2 



as expected from kinetic theory. 



Basic Statistical Mechanics 56 



Exercises 
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Boltzmann's postulate 



Recall that the Boltzmann distribution comes from 
Boltzmann's postulate 

S = fc B lnQ. 

The ability to derive the correct heat capacities for 
paramagnetic salts, ideal gases, and other systems are strong 
evidence that the postulate is correct. 

Results such as the heat capacity 

C = ^Nk B 

can be directly measured on real gases that are nearly ideal, 
such as Helium. 



This provides a means to obtain the Boltzmann constant, 
which turns out to be related to the gas constant: 

i. R 

N A 
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Some exercises 



Exercise 1 

Minimise the function f(x, y) = x 2 + y 2 , subject to the 
constraint x + y — 1 = 0. 

(i) Do it mentally, or by inspection of the graphs. 

(ii) Do it using Lagrange multiplier. 
Do the results agree? 
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Mentally ... 




1. The contours of f(x, y) are just circles, e.g. x 2 + y 2 = l, 

x 2 + y 2 = 2, ... The circles get larger as the value increases. 

2. The constraint x + y - 2 = 0 is a straight line with 
intercepts x = 2 or y = 2. 
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Using Lagrange multiplier ... 

1. Let g(x,y) = x + y - 2. The constraint is then g(x, y) = 0. 

2. The Lagrange function is 

f(x,y) + \g(x,y) = x 2 + y 2 + \{x + y - 2) (1) 

3. The gradient gives: 

1. Differentiate wrt x: 2x + \ = 0. 

2. Differentiate wrt y. 2y + A = 0. 

4. Solving with the constraint equation gives A = —2 and 
x = y = 1, as before. 
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Answer: Mentally 




1. By symmetry, f(x,y) is highest at the midpoint between the 
intercepts, where x = y = 1. 

2. Therefore, f(x, y) = l 2 + l 2 = 2. 
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Some exercises 



Exercise 2 

Use your calculator to work out In 10!. Compare your answer 
with the simple version of Stirling's theorem (iV In N — N). How 
big must N be for the simple version of Stirling's theorem to be 
correct to within (ii) 5% (ii) 1% ? 

Answer 

ln(10!) = 15.11 whereas 10ln(10)-10 = 13.03. There is less 
than 5% difference for N = 24 and less than 1% for B = 91. 
Stirling's approximation is quite accurate even for relatively 
small N. 
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Exercise 3 



Exercise 3 

Consider a 10000 distinguishable particles at room temperature, 
298 K. Suppose that each particle has 2 energy levels, 0.01 eV 
and 0.02 eV. Find the number of the particles in each energy 
level. ( Boltzmann's constant is 1.3807 x 10~ 23 J K -1 .) 
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Exercise 4 



Exercise 4 

Calculate the quantity k B T at room temperature, 298 K. Give 
the answer in eV. ( Boltzmann's constant is 1.3807 x 10~ 23 J 
K" 1 . Electron charge is 1.6 x 10" 19 C.) 

Answer 

In Joules, k B T = 1.3807 x 10" 23 x 298 = 4.114 x 10" 21 J. 

In eV, k B T/e = 4.114 x 10" 21 /l-6 x 10" 19 = 0.026 eV. 

Note: Since 1/40 = 0.025 is quite close to this answer, k B T at 
room temperature is often quoted as 1/40 eV. 
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Exercise 3 



Answer 

Since the number is quite large, we assume that the probability 
at each energy level is given by Boltzmann's distribution 

^exp(-e/fc B T). 

At 0.01 eV, Aexp(-s/k B T) = 0.6778A 
At 0.21 eV, Aexp(-s/k B T) = 0.4594A 
The total is 0.6778A + 0.4594.4 = 10000. So A = 8794. 

Therefore the numbers are: 

at O.OleV, 0.6778 x 8794 = 5961; 
at 0.02eV, 0.4594 x 8794 = 4040. 
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Exercise 5 



Exercise 5 

Consider a large number of distinguishable particles, at 
temperature T.. Each particle has 4 energy levels 0, k B T, 2k B T 
and 3k B T. Calculate the fraction of particles in each energy 
level. 

Answer 

Since the number is quite large, we assume that the probability 
at each energy level is given by Boltzmann's distribution 

Aexp(-e/fc B T). 
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First, calculate the following: 

at energy 0, exp(— e/k B T) = exp(O) = 1. 
at energy k B T, exp(-e/fc B T) = exp(-l) = 0.3679. 
at energy 2k B T, exp(-e/fc B T) = exp(-2) = 0.1353. 
at energy 3k B T, exp(-e/k B T) = exp(-3) = 0.0498. 

The sum is 1.553. The fractions are, therefore: 

at energy 0, 1/1.553 = 0.6439. 
at energy k B T, 0.3679/1.553 = 0.2369. 
at energy 2k B T, 0.1353/1.553 = 0.0871. 
at energy 3k B T, 0.0498/1.553 = 0.0321. 

Notice that fraction drops to nearly 1/10 after just 2k B T. 

In statistical mechanics, we often want to know if an energy 
level is likely to be populated at some temperature. If the level 
is much higher than k B T, then it is unlikely to be occupied. 
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Answer 



First, note that since Z = £ exp(-£i/k B T) 

N = Aj^expi-e./keT) = AZ 

i 

Next, 



d\nZ _ IdZ _ 1 „ e 
dT ~ ZdT ~ Z ^ ksT 



Multiplying the two equations above, 
8\nZ A 



N 



^== fcs ^E^exp(- £i /fc B T) 
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Exercise 6 

Given that the total energy of a system is 

U = AJ2^ exp(-e z /k B T) 

i 

where A is defined by the normalisation from 

N = Aj^expi-eJkBT) 

i 

demonstrate that 

9 <9lnZ 
U = Nk B T 

where 

Z = ]rexp(- £i /fc B T) 

i 
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Rearranging, 



Nk B T 2d ^- = Aj2eiexp(-sJk B T) 



which is U. Therefore 



7 d\nZ 

Nk B T = U 
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Exercise 7 

The total energy of a paramagnet is 

U = -N)i B Btanh(n B B/k B T) 

Show that the heat capacity is given by 

\ 2 exp(0/T) 



C = Nk B 
where 0 = 2^ B B/k B 



TJ (exp(0/T) + l) 2 



Show that the limiting form at high temperature is 

and at low temperature is 

C = Nk B Q 2 exp(-0/T) 
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At high temperature, exp(0/T) — > 1. 

?6\ 2 exp(0/T) 
vtJ (exp(0/T) + l) 2 

W x2 1 



C = iVfc B 



TJ (l + l) 2 



iVfc B /0^ 2 



4 VT7 

At low temperature, exp(0/T) + 1 -> exp(0/T). 

2 exp(0/T) 



C = iVfc B 



77 (exp(0/T) + l) 2 
0\ 2 exp(0/T) 
T7 exp(0/T) 2 

= JVfcs(^) exp(-0/T) 
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Answer 



The heat capacity is obtained by differentiating 

U = -Nfj, B Btanh(n B B/k B T): 



c = % = - N ^ B \MTi) 



cosh (>**) 



In terms of 6 = 2^ B B/k B , 



C = Nk R - 



\2Tj 



-V 



cosh 



= Nk B l — ) 



= Nki 



= Nk B 



(IX 4 

V2T7 (exp(0/2T) + exp(-0/2T)) 2 
0\ 2 l exp(0/T) 

TJ (exp(0/2T) + exp(-0/2T))2'exp(0/T) 
2 exp(0/T) 



77 (exp(0/T) + l) 2 
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Exercise 8 

Show that the partition function of an ideal gas: 

/-co 4-mivV i /o 

Z= / ^-(2m £ ) 1 / 2 exp(-e/lt B T) ( i£ 

is given by: 



Z = V 



/2irmk B T\ 3 / 2 
I h* ) ■ 



[You are given that 

° 1 ,2n (J -a 2 x 2 ^ T _ 

n!(2a) 2 »a 



J— CO 
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First, move the constant factors outside the integral: 

Z = 4m f(2m) 1 /2 r e l/2exp(- e /fc B T)d e 
h 3 JO 

There is no power of 1/2 in the integration formula that is 
given. So we start with the substitution 



e = x 2 



Then 



Z = =— (2m) 1 ' 2 / xexp(-x z /k B T).2xdx 



h 3 Jo 

r ( 2lll )' 2 2 J™x 2 exp(-x 2 /k B T)da 



h 3 

Next, set n to 1 in the given formula 



» i 2 2 27T 1 / 2 

oo (2a) z a 



The lower of limit of this integral is — oo. To make this agree 
with the integral for the partition function Z ... 
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We get 



Z = 4 ^(2 TO )V2 2 !^ (fcBT) 3/2 



Simplifying, 

'2mnk B T\ 3 / 2 
1? J 

we get the formula for the partition function. 



Z = V 
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We change the lower limit to 0 and halve the right hand side, 
and the formula becomes: 

(2a) 2 a 

Now comparing with the partition function again: 

ArmvV xi/o_ f°° 



roc 2 2 

/ x 2 e~ a x dx 
JO 



Z = 



-(2m) 1 / 2 2 J°°x 2 exp(-x 2 /k B T)dx 



h 3 

the integrals become the same if we set 



a 2 = 



k B T 

Substituting this into the formula gives 



r x 2 e -* 2 /k B T dx = ^ kBT) 3/2 
JO 4 



Substituting this into partition function intergral, 
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Exercise 9 

What is the root-mean-square speed of a helium atom at room 
temperature? 

[Given that the relative atomic mass of helium is 4. 
Atomic mass unit u is 1.6605 x 10~ 27 kg. 
Boltzmann's constant k B is 1.3807 x 10~ 23 J K -1 . ] 
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First, we find the mean of v 2 . This is can be obtained by 

- 1 f°° 2 

v 2 = — / v n(e)de. 
NJo w 

We have multiplied v 2 by the number particles at each energy 
interval de, and divided by the total number N. 

The energy of the atom in an ideal gas is just the kinetic 
energy, 

1 2 

e = —mv . 

2 

If we rewrite the above integral for v 2 in terms of the energy, 
we get 

7} ^ /*°° 

v 2 = / en(e)de. 

This integral is just the total energy U, i.e. 

v 2 = —U. 
Nm 
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We can now calculate the root-mean-square speed. The mass 
of the helium atom is 

Au = 4 x 1.6605 x 10~ 27 = 6.642 x 10~ 27 kg. 
Taking the room temperature as 298 K, the answer is 



3k B T 



3 x 1.3807 x 10" 23 x 298 
~ \ 6.642 x 10" 27 

= 1363 m/s 
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We have previously obtained the energy U from the partition 
function: 

U = ^Nk B T. 
So the mean square speed is 

-o 2 3„„ _ 3k B T 
v 2 = Nk B T = — — 

The root-mean-square speed is then 

3k B T 
m 
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